Stability of coreless vortices in ferromagnetic spinor Bose-Einstein condensates 
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We study the energetic and dynamic stability of coreless vortices in nonrotated spin-1 Bose- 
Einstein condensates, trapped with a three-dimensional optical potential and a Ioffe-Pritchard field. 
The stability of stationary vortex states is investigated by solving the corresponding Bogoliubov 
equations. We show that the quasiparticle excitations corresponding to axisymmetric stationary 
states can be taken to be eigenstates of angular momentum in the axial direction. Our results show 
that coreless vortex states can occur as local or global minima of the condensate energy or become 
energetically or dynamically unstable depending on the parameters of the Ioffe-Pritchard field. The 
experimentally most relevant coreless vortex state containing a doubly quantized vortex in one of 
the hyperfine spin components turned out to have very non-trivial stability regions, and especially 
a quasiperiodic dynamic instability region which corresponds to splitting of the doubly quantized 
vortex. 

PACS numbers: 03.75.Mn, 03.75.Kk, 67.57.Fg 



I. INTRODUCTION 



Development of optical trapping techniques for alkali 
atoms has enabled experimental studies of dilute atomic 
Bose-Einstein condensates (BECs) with spin degrees of 
freedom [j], For these systems, the order parame- 
ter is a spinor field which can exhibit a rich variety of 
different topological textures ranging from coreless vor- 
tices @, H| to monopoles [f| @ and Skyrmion-type con- 
figurations @, H, Q. In scalar BECs, a vortex is always 
fully characterized by the phase winding about the vor- 
tex core, whereas in spinor BECs the characteristics of 
a vortex are determined by winding numbers of different 
components as well as the core polarization - for coreless 
vortices the superfluid velocity is non-divergent at the 
vortex core and hence the vortex core can be polarized. 

Coreless vortices are topologically unstable, i.e., they 
can be continuously deformed to a uniform texture. 
Thus their existence as stable states typically requires 
the presence of additional forces such as interactions at 
large distances from the vortex core or external fields 
which impose nontrivial asymptotic boundary conditions 
fiol . [ri| . Coreless vortices such as the Mcrmin-Ho vor- 
tices in spinor BECs are analogous to those in superfluid 
3 Hc-A [3J, in which they appear as equilibrium objects 
if the system is rotated externally. Thus it is natural 
to assume that such objects would be generated also in 
rotated gaseous condensates. It has indeed been theoreti- 
cally confirmed that also in these systems the Mcrmin-Ho 
vortices are energetically stable for certain values of the 
trap rotation frequency and magnetization 

Manipulating spinor condensates with external mag- 
netic fields has been in vogue among both theorists and 
exp erimentalist during the recent years [l3|, [Hi HE 
Il7| . Topological phase engineering by time-dependent 
external magnetic fields has been successfully used to 



create vortex structures 

[HHSl. Recently, 
Leanhardt et al. succeeded in creating a coreless vortex 
in a F = 1 spinor condensate in a Ioffe-Pritchard (IP) 
magnetic trap [22| by adiabatically switching off the mag- 
netic bias field along the trap axis [l5| ■ The ground state 
phase diagram corresponding to the IP field combined 
with an optical confinement potential has been computed 
for a condensate uniform in the direction of the vortex 
axis [l|| , and it shows that the IP field renders the ground 
state of the system to be a coreless vortex. This is due to 
the tendency of the spin to align with the external field, 
which leads to formation of a coreless vortex. An example 
of the spin texture in an IP field is shown in Fig. Q] in the 
pancake-shaped geometry of Ref . [l6| . Zhang et al. [l~7| 
found that within an adiabatic approximation (l3j |. the 
difference of atomic spatial angular momentum and hy- 
perfine spin is conserved in the IP field. This conservation 
law implies that the ground state of the condensate car- 
ries a persistent current with definite winding numbers. 




FIG. 1: (Color online) One possible spin texture of a coreless 
vortex in pancake-shaped condensate in the presence of the 
Ioffe-Pritchard field. Spin of the condensate tends to align 
with the external field. 

In the previous studies [H,[l3| only the global energetic 
stability of coreless vortex states was considered. How- 
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ever, the coreless vortices created in the experiment |15| 
contain a doubly quantized vortex in one of the hyper- 
fine spin states, and it is of interest to find out whether 
such coreless vortex states inherit the dynamic instabil- 
ity of doubly quantized vortices in scalar condensates 
IS HI M |H HI- Furthermore, the topological 
phase imprinting methods [TJ, EH, US H3, HH used to 
create doubly quantized vortices involve adiabatic inver- 
sion of the bias field in the presence of the IP field. The 
intermediate state in this process is always a coreless vor- 
tex which for certain field parameters is the ground state 
of the system but otherwise its stability properties are 
so far unknown. Hence the study of stability of coreless 
vortices can also give valuable insight to the process of 
vortex creation using topological phase engineering. 

In this paper, we analyze further the vortex stabil- 
ity phase diagram of ferromagnetic spinor condensates in 
the presence of a IP field and an optical trap potential 
by computing the quasiparticle spectra from the Bogoli- 
ubov equations in finite three-dimensional trap geome- 
tries, for oblate, spherical, and prolate condensates. We 
show that different coreless vortex states are locally ener- 
getically stable for a wide range of IP field configurations, 
although for other parameter values dynamic and ener- 
getic instabilities may occur. Quasiparticle states corre- 
sponding to these instabilities reveal that the dynamic 
instability of coreless vortices containing doubly quan- 
tized vortices in one of the hyperfine spin components 
is similar to the dynamic instability of the doubly quan- 
tized vortices in scalar condensates, that is, the doubly 
quantized vortices tend to split due to the instability. 

II. MEAN FIELD THEORY 



Above, m is the mass of the atoms, \x the chemical po- 
tential, and Co, c 2 are the coupling constants related to 
s-wave scattering lengths in different total hyperfine spin 
channels Q . Depending on whether the interaction cou- 
pling constant C2 is positive or negative, the condensate 
is either antiferromagnetic or ferromagnetic, respectively. 
Klausen et al. [2i| have shown that 87 Rb is ferromag- 
netic, whereas 23 Na atoms realize an antiferromagnetic 
condensate [3(|. We take the confining potential U(r) 
created by an external optical field to be the axisymmet- 
ric harmonic potential 

where A = u> z /u) r . 

In the presence of an additional external magnetic field 
B(r), we take also into account the linear Zeeman term 
/i B 3FB(r) • S(r) in the energy functional. The constant 
fi B is the Bohr magneton and g F is the Lande g-factor. 
The free energy becomes in this case 

FbP) = S[i>] +Mb5f / drB(r) ■ S(r), (3) 

and the magnetic field is in our analysis of the Ioffc- 
Pritchard form 

B(r) = B ± (xx-yy)+B z z. 

Stationary states of the condensate satisfy 8T[ip]/ 'Sip = 0, 
which yields the Gross-Pitaevskii (GP) equation 



In the mean-field approximation, we describe the con- 
densate formed by weakly interacting ultra-cold bosonic 
atoms in the z-quantized basis \F = l,m F ), m F = 
— 1,0, 1, with a spinor field ip = (^1^0 ip-i) T and the 
free energy functional of the form y 



dr 
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|V^(r)| 2 + C/(r)|V(r)| 2 -/i|V(r)p 
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where the spin density is given by 

S(r) = 5>:(r)F ai ,V6(r) 



(1) 



(2) 
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and the angular momentum matrices F = [F x F y F Z ) T 
are the usual generators of the spin rotation group 5*0(3) 
with matrix representations 



H[i>]i> = fiifi, (4) 
where the non- linear operator Ti.[ijj] is given by 

2rr 

c 2 F • S(r). 



n[ip] = - -^V 2 + U(r) + ^g F B(r) ■ F + c \^(r)\ 2 
Zm 



The quasiparticle spectrum corresponding to a given 
stationary state can be solved from the generalized Bo- 
goliubov equations [ll, HH 
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v q (r) 



flUJn 



u q (r) 
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where u q = {u ql u q0 u^f and v q = {v ql v qfi V g _i) T 
are the quasiparticle amplitudes and the operator T> is 
defined as 
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such that the components of the A and B operators are 
Aij =(- |^V 2 +U(t)- f^jSij + /z B < ?P ^ J B a (r)(F a ) ij 



a k, I 

Bij =coipi(r)i/)j(r) + c 2 ^2^2{F a ) lk (F a ) J i^ k (r)^i(r), 

a k, I 

for i,j <= {1,0,-1}. Due to the conjugate symmetry 
of the Bogoliubov equations, we may concentrate only 
on the quasiparticle modes for which the quadratic form 
/ dr(|u g (r)| 2 — |w g (r)| 2 ) is non-negative. The quasipar- 
ticle spectrum determines the stability properties of the 
corresponding stationary state: If the quasiparticle spec- 
trum contains excitations with non-real frequencies u q , 
the system is dynamically unstable and even small initial 
perturbations can render the population of these modes 
to grow exponentially in time even in the absence of dis- 
sipation. On the other hand, if there exists modes with 
negative eigenfrequencies, the state is energetically un- 
stable, and in the presence of dissipation the condensate 
can lower its energy by transferring particles from the 
condensate state to such anomalous modes. Vice versa, 
if all the quasiparticle eigenfrequencies are positive, the 
state is locally energetically and dynamically stable, and 
should be long-living and robust against small perturba- 
tions even in the presence of dissipational mechanisms. 



III. AXISYMMETRIC VORTEX STATES 

Within an adiabatic approximation, that is, assuming 
that hyperfme spins of the atoms remain aligned (or an- 
tialigned) with the local magnetic field, the ground state 
of the condensate in the IP field is axisymmetric p7| . 
However, nonadiabatic effects are important to some ex- 
tent, since in Rcf. fl(| the ground state phase diagram 
was shown to contain a non-axisymmetric vortex state 
in the antiferromagnetic case. On the other hand, we 
have verified with fully three-dimensional (3D) compu- 
tations, i.e. without setting any symmetry restrictions to 
the condensate wavefunction, that for ferromagnetic con- 
densates with nonzero radial IP field strength, the ground 
state is always axisymmetric for all the parameter values 
considered. Consequently, we restrict to consider only 
configurations that are axisymmetric. 

In cylindrical coordinates r = (r, ip, z), axially sym- 
metric vortex states are of the form 

^P(v) = (%l>i(r,z)e lK ^ iP (r 7 z)e lKoip V-i (r, z)e iK -^) T , 

(6) 

where Ki G Z are the winding numbers of the 
three components — we refer to states of this form as 



Ko, K-i). In addition, all physically measurable den- 
sities corresponding to axisymmetric states have to be 
axisymmetric. When applied to spin-density, this re- 
quirement combined with Eqs. @ and ^ implies the 
relation 



2kq = Ki + K—X, 



(7) 



between the winding numbers. The angular dependence 
of the IP field sets additional restrictions for the ground 
state to be axisymmetric: For the Zeeman energy term 
/! B g F B(r) • S(r) to be rotationally symmetric the wind- 
ing numbers of the condensate state have to satisfy the 
additional constraints 



K() = Kl — 1 



1. 



(8) 



One notes that these latter relations imply also the rela- 
tion in Eq. ([7]), and are thus more restrictive. We have 
verified with fully 3D computations without any sym- 
metry assumptions that the ground state of system in- 
deed satisfies the restrictions given in the Eqs. 
The states with the lowest angular momenta satisfying 
Eqs. (7J) and © are (2,1,0), (1,0,-1), and (0,-1,-2) 
which are coreless vortex states. The numerical results 
indeed show that the ground state in the presence of a 
strong enough IP field contains a coreless vortex. 

The coreless vortices in the states (2,1,0) and 
(0,-1,-2) are ferromagnetic in the sense that the spin 
of the condensate is aligned (or antialigned) with the ex- 
ternal field also in the core region. Furthermore, (2, 1, 0) 
and (0, —1, —2) are equivalent since they differ only by in- 
version of the spin quantization axis. On the other hand, 
(1, 0, —1) is polar in the sense that at the vortex core S(r) 
vanishes [1, [3l|. The difference in the spin texture turns 
out to be significant for the phase diagram and the local 
energetic stability of the vortex state. Due to vanishing 
spin density at the vortex core, the spin texture of the 
polar vortex (1, 0, —1) can align with the IP field equally 
well for both positive and negative values of the bias field 
B z . 

In investigating the stability properties of stationary 
states by solving the Bogoliubov equations, it is to be 
noted that apart from the possible degeneracy of spec- 
trum, axisymmetric states can have non-axisymmetric 
quasiparticle excitations also due to the fact that T> does 
not in general commute with the angular momentum op- 
erator L z = ~ihd v . We can, however, show that if the 
winding numbers of different components satisfy Eq. l[8]). 
then there exists a unitary transformation U such that 
the transformed Bogoliubov operator T> = W~DU com- 
mutes with L z . The unitary transformation IA is in this 
case given by 



U = 





u\ 



where U Q = 



e' KlV 
e iK °f 











for the Bogoliubov operator V = T>[t/j] corresponding to 
an axisymmetric state {k\, kq, Now [T>,L z ] = 
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follows from the condition stated in Eq. JH]). Hence the 
eigenstates of V can be chosen to be eigenstates of L z . 
Writing the eigenvalue equation T>w = rjw in the form 
UVU^Uw = rjUw and taking into account that UVU^ = 
X>, we observe that the quasiparticle amplitudes in the 
original Bogoliubov equation (O can be taken to be of 
the form 

u q ^v)=u^(r,z)e^ +K >, (9) 
v q j(r) = v qJ (r,z)e i ^- K ^, i = -1,0,1, (10) 

where K q is an angular momentum quantum number of 
the excitation. In addition to this analytical argument, 
we have verified numerically without any symmetry as- 
sumptions that for the pancake-shaped condensates with 
A = uj z /uj r 1, all the low-energy quasiparticle states 
are axisymmetric. 



IV. NUMERICAL RESULTS 

In the following we consider ferromagnetic condensates 
which can be realized, e.g., with 87 Rb atoms. For the 
scattering lengths of 87 Rb, the estimate of van Kempcn 

-0.005. In this 



et al. 32 1 implies the ratio c^jc^ 



study we use the value c 2 /c = —0.02 which was also 
used in the previous theoretical studies of Mizushima et 
al. [121 [33( 1 . In the numerical calculation we use dimen- 
sionlcss units in which cq = cqtuN / a r h, where the charac- 
teristic length scale of the trap is given by a r = \Jhjmio r . 
We adopt cq — 10000, which corresponds, for example, 
a condensate of N = 8.5 x 10 5 87 Rb atoms trapped 
such that u! r = 2ir x 200 Hz. We take the Lande g- 
factor to be g F = — ^ which indicates that the hyper- 
fine spine tends to align with an external field. We have 
searched for solutions of the GP equation and the Bo- 
goliubov equations by using finite-difference discretiza- 
tion combined with relaxation methods and the implicitly 
restarted Arnoldi method implemented in the ARPACK 
numerical library [34j]. In the numerical calculations, we 
also use the conjugate symmetry of the Bogoliubov equa- 
tions and consider only positive K q , which slightly reduces 
the numerical effort. 

Numerical computations showed that the ground state 
configuration for non-vanishing perpendicular IP field 
B± is always one of the coreless vortex states (2,1,0), 
(1,0,-1) or (0,-1,-2). Figure [2] shows the computed 
stability phase diagrams displaying the dynamic and en- 
ergetic stability /instability regions for the state (2,1,0) 
as a function of the IP field parameters. The correspond- 
ing phase diagram for the state (1,0,-1) is shown in 
Fig. El The results have been computed for the trap 
asymmetry parameter values A = 0.2, 1.0, and 5.0, cor- 
responding to prolate, spherical, and oblate geometries, 
respectively The phase diagram for the state (0, —1, —2) 
is the same as for (2, 1, 0) if the sign of the bias field B z 
is reversed. 
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FIG. 2: Stability phase diagram for the axially symmetric 
(2,1,0) vortex state. On the top A = 0.2, in the middle 
A = 1.0, and at the bottom A = 5.0. The left panels show the 
different stability phases and the right panels the maximum 
value of |Im{k; 9 /(Av}| for each point in the (B±,B Z ) plane 
as a grayscale plot: bright regions correspond to dynamically 
stable regions, and dark ones dynamically unstable. In the re- 
gion denoted as "local minimum" (2, 1,0) is a local minimum 
of the mean field energy and the polar state (1, 0, —1) is the 
ground state of the system. The dynamic instability region 
indicates the section of the parameter space where dynamic 
instability modes can appear. 




FIG. 3: Stability phase diagram of the polar vortex (1, 0, —1) 
(right panel) and the lowest energy excitations as a function 
of B z for fixed jj, s B±/2huj r — 0.5 (left panel). The phase 
diagram is symmetric with respect to reversing the sign of 
B,. 



From the diagrams one infers that the ferromagnetic 
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FIG. 4: Quasiperiodic structure of the dynamic instability as a function of B z for the axially symmetric (2, 1,0} vortex and 
= 0.7. Left panel: A = 5.0; middle: A = 1.0; and right panel: A = 0.2. 



state (2,1,0) and the polar state (1,0,-1) can co-exist 
as local energy minima for suitable IP field parame- 
ters. The transition between the ferromagnetic and polar 
ground states is determined by the relative magnitude of 
the IP field components: for B± ^> \B Z \ the ferromag- 
netic core of the vortex states (2,1,0) and (0,-1,-2) 
becomes energetically unfavorable, rendering the polar 
state (1,0,-1) the global minimum of energy, and vice 
versa for \B Z \ 3> B±. 

The negative energy anomalous modes indicating en- 
ergetic instability occur for quantum numbers K q = ±1 
for the polar state (1, 0, —1) and for n q = ±1, ±2 for the 
ferromagnetic state (2,1,0). The (1,0,-1) state turned 
out to be dynamically stable for all the parameter values 
investigated, but the states (2, 1,0) and (0, —1, —2) have 
complicated dynamic instability regions. The dynamic 
instability modes occur only for k q = ±2. Figure [2] shows 
the greyscale plots of the maximum imaginary parts of 
the Bogoliubov eigenfrequencies for the state (2, 1,0) as 
functions of the IP field parameters for three trap geome- 
tries. We note that the regions marked as dynamically 
unstable in Fig. [2] contain narrow stripe-like patterns in 
which the vortex state is in fact dynamically stable. How- 
ever, especially for the prolate geometry these stripes are 
very narrow and probably experimentally indistinguish- 
able. The stripe-like quasiperiodic structure of the dy- 
namic instability can be observed more clearly in Fig. 2J 
in which the maximum imaginary part is plotted for a 
fixed value of B±. One observes that the magnitude of 
the largest imaginary part oscillates markedly before it 
saturates for strong enough bias fields B z . 

Figure [5] shows that the (2,1,0) vortex tends to be- 
come dynamically more stable with increasing A, i.e., the 
vortex is generally more stable in the pancake shaped 
geometry than in the cigar shaped one. This is due 
to the suppression of the density of low-energy excita- 
tions in the limit of tight confinement in the z direc- 
tion. However, the dynamic instability persists even for 
A = uj z /uj 7 . ^> 1. The energetic stability of the coreless 
vortices is intimately related to dynamic instability since 
it has been shown in Ref. [28[ that dynamic instabilities 
in scalar condensates are formed when a negative and a 



positive energy excitation are in resonance. Similar phe- 
nomenon takes place also in spinor condensates, and can 
be observed in this case by inspecting the spectrum of 
the quasiparticle energies for n q = ±2 shown in Fig. [5] 
For A = 1.0 and fi B B±/2huj r = 0.5, the dynamic instabil- 
ity may occur for parameter values fi B B z /2huj r > 1.16. 
From Fig. [5] one observes that some positive energy ex- 
citations are missing in the region [i B B z /2huj r > 1.16 
due to the resonance. We note that the regions in which 
the (2,1,0) vortex is dynamically stable seem to be man- 
ifested by the absence of certain positive and negative 
energy excitations. For the polar vortex (1,0,-1) the 
excitation spectrum is much simpler containing only one 
anomalous mode. The stability phase diagram for the 
polar state (1,0,-1) in spherical geometry is shown in 
Fig. E] together with the lowest quasiparticle energies. 
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FIG. 5: The negative energy excitations and the correspond- 
ing positive energy excitations with K q = ±2 for the ferro- 
magnetic vortex state (2,1,0). Here HuB±/2ftjjj r — 0.5 and 
A = 1.0 are fixed. The dynamic instabilities start to occur 
for [i-bB z /2hu) r > 1.16 which is denoted by the solid line. 
The intervals denoted by the dashed lines and the Roman 
numerals I, II, III, and IV correspond to regions of values of 
B, for which the (2,1,0) vortex state is dynamically stable 
(cf. Fig.[2j. 

For the ferromagnetic vortex state (2,1,0), the unit 
vector field n(r) = S(r)/|S(r)| is well-defined everywhere 
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inside the cloud. Thus the spin texture of this state can 
be mapped to a simply connected subset S C S 2 and the 
perpendicular part of the IP field prevents the subset S 
from shrinking to a point. The surface S' corresponding 
to the spin texture of the polar vortex (1, 0, —1) contains 
a hole due to vanishing S(r) at the vortex core. Thus 
the two spin textures arc topologically incquivalent and 
there is a topological phase transition of the ground state 
between the two textures occurring at the boundary be- 
tween the regions of global and local stability in Fig. [51 
To investigate qualitatively the nature of the dynamic 
instability, we consider states of the form 

$(r) =V(r)+7?K(r)+u*(r)], (11) 

where a slightly populated quasiparticle state corre- 
sponding to the dynamic instability with the largest 
imaginary frequency has been added to the correspond- 
ing stationary state. Isosurfaces of particle densities in 
different hypcrfinc spin components of such slightly per- 
turbed (2, 1, 0) vortex states are shown in Fig. [5] In the 
component containing the doubly quantized vortex, one 
observes a helical vortex chain structure, which is very 
similar to the one discovered in the numerical calcula- 
tions of Huhtamaki et al. [25| for the splitting times of 
doubly quantized vortices in scalar condensates. Thus 
we observe that for the (2, 1,0) vortex, the dynamic in- 
stability is essentially due to the splitting of the doubly 
quantized vortex in the m F — +1 component. Due to the 
helical structure of the splitted vortex in this multicom- 
ponent case, the splitting of the doubly quantized vortex 
is accompanied with a spiral flow about the z-axis of the 
vortex free component m F = — 1. 




m F = — 1 m F — m F — 1 



FIG. 6: (Color online) Isosurfaces of particle densities in dif- 
ferent spin components corresponding to vortex state (2, 1,0) 
and a slight excitation of a dynamic instability mode. The 
trap asymmetry parameter is A = 0.2. Population of the 
dynamic instability mode is 1% of the total number of par- 
ticles. The isosurfaces correspond to values | SP i (r) ] 2 = 
2.6 ■ PT 5 N/al, |*o(r)| 2 = 1.0 ■ KT 5 N/a z r , and |*-i(r)| 2 = 
1.0 • 10~ 5 N/al. 

Due to the Zeeman energy term, the average spin S(r) 
tends to align with the local magnetic field. However, in 



the vicinity of the vortex core the average spin may devi- 
ate from the direction of the local magnetic field without 
costing too much energy. The spin texture of the state 
(2, 1,0) outside the vortex core points always along the 
local magnetic field, but for B z > 0, S(r) has opposite 
direction to the magnetic field. Since the perpendicular 
part of the IP field is topologically nontrivial with wind- 
ing number — 1 [15l . |35| , a continuous transformation of 
the spin texture to energetically more favorable state can- 
not be accomplished without spin rotations against the 
local magnetic field. Thus the IP field creates an energy 
barrier which prevents the spin texture from unwinding 
itself via spin rotations. 

V. CONCLUSIONS 

In conclusion, we have studied the local stability of 
corcless vortex states in optically trapped ferromagnetic 
spinor F = 1 BECs in the presence of the Ioffc-Pritchard 
field. The ground states of the system turned out to 
be axisymmetric and it was shown analytically that the 
corresponding quasiparticle states are also axisymmet- 
ric. By computing numerically the quasiparticle spectra 
we have shown that there can co-exist different coreless 
vortex states as local and global minima of energy for 
a wide variety of different external field configurations. 
The experimentally most tractable vortex configuration, 
the coreless vortex state (2, 1, 0), was shown to possess a 
rich phase diagram in which the vortex transforms grad- 
ually from a global minimum of energy to a dynamically 
unstable stationary state as the bias field of the Ioffe- 
Pritchard trap is ramped up from negative to positive 
bias. Based on these results one should be able to real- 
ize experimentally robust coreless vortex states by load- 
ing the atoms into an optical trap with an IP field and 
allowing the condensate to relax into the ground state. 
This method should also enable experimental creation of 
the polar vortex state (1,0,-1). Interesting questions 
for the future research would be to investigate the phase 
diagram of the ground state and the local stability of 
vortex states in an external Fricdburg-Paul (hexapole) 
magnetic field (36l . \3l\ , and the exact decay mechanisms 
of energetically unfavorable spin textures. 
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